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Abstract 

We consider nondecreasing entropy solutions to 1-d scalar conservation laws and 
show that the spatial derivatives of such solutions satisfy a contraction property 
with respect to the Wasserstein distance of any order. This result extends the L 1 - 
contraction property shown by Kruzkov. 



Existence and uniqueness of solutions to scalar conservation laws in one space dimen- 
sion have been established by Kruzkov in the framework of entropy solutions (see |3] for 
instance), and among the properties satisfied by these solutions it is known that the L 1 
norm between any two of them is a non-increasing function of time. 

In this work we shall focus on a class of entropy solutions such that a certain distance 
between the space derivatives of any two such solutions is also nonincreasing in time. On 
this class of solutions this result extends the L 1 norm contraction property. 

More precisely we consider as initial data nondecreasing functions on M. with limits 
and 1 at — oo and +00 respectively. These properties are preserved by the conservation 
law, and corresponding solutions have been shown in [2j to arise in some models of pres- 
sureless gases, obtained as a continuous limit of systems of sticky particles. Noticing that 
the distributional space derivative of these functions are probability measures, we may 
consider the Wasserstein distance between the space derivatives of any two such solutions, 
and we shall prove in this paper that this distance is a nonincreasing function of time, 
constant in the case of classical solutions. 
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1 Introduction to the results 

Given a locally Lipschitz real-valued function / on R, called a flux, we consider the 
scalar conservation law 



(1) 



u t + f(u) x = 0, t > o, x e 
w(0,.) = u°, 



with unknown u = u(t,x) G R and initial datum tt° G L°°(R), and where the subscripts 
stand for derivation. 

We shall consider solutions that are called entropy solutions (see [5J for instance) and 
are defined as follows: a function u = u(t,x) G L°°([0, +oo[xR) is said to be an entropy 
solution of (II} on [0, +oo[xR if the entropy inequality 

(2) E{u) t + F(u) x < 

holds in the sense of distributions for all convex Lipschitz function E on R, and with 
associated flux F defined by 

(3) F(u) = P f(v)E'(v)dv. 

Jo 

This means that 

(E(u) (ft + F(u) <p t ) dtdx + / E(u°(x)) ip(Q, x) dx > 

Jr 

for all nonnegative (p in the space C£°([0, +oo[xR) of C°° functions on [0,+cxd[xR with 
compact support. 

We shall also consider classical solutions, that is, functions u = u(t, x) in C 1 (]0, +oo[xR)n 
C([0, +oo[xR) satisfying ^ pointwise. 

In particular any classical solution to Jl} satisfies (J2J), i.e. is an entropy solution, and 
conversely any entropy solution satisfies (JH) in the distribution sense. 

For entropy solutions, the following result is due to Kruzkov (see 

Theorem 1.1 For every u° G L°°(R) ; there exists a unique entropy solution u to (J2J) in 
L°°([0, +oo[xM) n C([0, +oo[, L} oc (M)). 

Moreover for classical solutions, we have (see [H] for instance): 

Theorem 1.2 Given a C 2 flux f and a C 1 bounded initial datum u° such that f o u° is 
nondecreasing on R, the unique entropy solution u to Q is a classical solution. 

In this work we shall consider initial data in the subset U of L°°(R) defined by 
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Definition 1.3 A function v : R — > R belongs toll if it is nondecreasing, right- continuous, 
and has limits and 1 at — oo and +oo respectively. 

The following proposition expresses that this set is preserved by the conservation law 

®: 

Proposition 1.4 Given an initial datum u° 6 £/ie entropy solution u given by Theorem 
li.il «s sac/i t/iat n(t, .) belongs to U for all t > 0. 

More precisely, given any t > 0, the L°°(R) function n(£, .) is a.e. equal to an element 
of the set U, which on the other hand is characterized by 

Proposition 1.5 The distributional derivative v x of any v G U is a Borel probability 
measure on R, and for any 

v(x) = v x (] - oo,x]). 
Conversely, if fi is a probability measure on R, then v defined on R as 

v(x) = [/,(] — oo,x]) 

belongs to U, and v x = fi. 

Consequently the map v h- > v x is one-to-one from U onto the set V of probability 
measures on R (and U can be seen as the set of repartition functions of real-valued 
random variables). 

Propositions 11.41 and 11.51 allow us to characterize at any time the distance between 
two solutions (with initial datum in U) in terms of their space derivatives, in particular 
by means of the Wasserstein distances: given any real number p > 1, the Wasserstein 
distance of order p is defined on the set of probability measures on R by 

W p (fi, p) = inf ( [ \x- y\ p dn(x, y) \ 

where 7r runs over the set of probability measures on R 2 with marginals \i and jl; these 
distances are considered here in a broad sense with possibly infinite values. 

This paper aims at proving that the Wasserstein distances between the space deriva- 
tives of any two such entropy solutions is a nonincreasing function of time: 

Theorem 1.6 Given a locally Lipschitz real-valued function f on R and two initial data 
u° and u° in U, let u and u be the associated entropy solutions to JH). Then, for any t > 
and p > 1, we have (with possibly infinite values) 

W p (u x (t,.),u x (t,.)) < W p (ulu° x ). 
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We shall see in Section [2] that for p = 1 the distance W\ satisfies 

Wi{v x ,v x ) = \\v - V\\ L 1(M) 
for all v, v G U. Hence Theorem II .61 reads in the case p — 1: 

\\u{t,.) - u(t,.) ||ii( R ) < ||U°- W°||l,i(R). 

Thus, for initial profiles in U, we recover the /^-contraction property given by Kruzkov. 

In the case of classical solutions, the result of Theorem 11.61 is improved, since the 
Wasserstein distance between two solutions is conserved: 

Theorem 1.7 Given a C 1 real-valued function f on M., let u° and u° in U be two initial 
data such that the associated entropy solutions u and u to JH) are classical solutions, 
increasing in x for allt > 0. Then for any t > and p > 1 we have (with possibly infinite 
values) 

W p (u x (t,.),u x (t,.)) = W p (u° x ,u° x ). 

>From these general results can be induced some corollaries in the case of initial data 
in the subsets U p of U defined as: 

Definition 1.8 Let p > 1. A function v inlA belongs tolA p if its distributional derivative 

v x has finite moment of order p, that is, if I \x\ p dv x (x) is finite. 

Jr 

As in Proposition 11.51 the map v i— > v x is one-to-one from U p onto the set V p of 
probability measures on M with finite moment of order p. But we shall note in Section 
El that the map W p on V p x V p defines a distance on V p . Then the real-valued map d p 
defined on U p x U p by 

d p (v,v) = W p (v x ,v x ) 
induces a distance on U p , and for the associated topology we have 

Corollary 1.9 Given a locally Lips chitz function f on E, p > 1 and u° G U p , the entropy 
solution u to (JIJ belongs to C([0, +oo[, IA P ). 

In particular for p — 1 

di{v,v) = Wi(v x ,v x ) = \\v - v\\ L i(R), 
and the previous result can be precised by 

Corollary 1.10 Given a locally Lipschitz function f on R and u° G U\, the entropy 
solution u to ((IJ is such that 

\\u(t, .) - u(s, .)||li( R ) < \t ~ s\ ||/'||l°°([0,1])- 
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This known result holds under weaker assumptions (for u° with bounded variation, see 
but in our case it will be recovered in a straightforward way. 

Finally Theorem 11.71 can be precised in the U p framework in the following way: 

Corollary 1.11 Given a C 2 convex flux f and two C 1 increasing initial data u° and u° 
in Up for some p > 1, the following three properties hold: 

1. the associated entropy solutions u and u are classical solutions; 

2. u(t, .) and u(t, .) belong to U p and are increasing for all t > 0; 

3. for all t > 0, we have (with finite values) 



The paper is organized as follows. The definition and some properties of Wasserstein 
distances are discussed in greater detail in Section [2j In Section 03 we consider the case 
of classical solutions, proving Theorem 11.71 and Corollary II .111 Then the general case of 
entropy solutions is studied in Sections |U and El more precisely in Section HI we introduce 
a time-discretized scheme, show the W p contraction property for this discretized evolution 
and prove the convergence of the corresponding approximate solution toward the entropy 
solution; Theorem 11.61 and its corollaries follow from this in Section In Section M we 
shall finally see how such results extend to viscous conservation laws. 

2 Wasserstein distances 

In this section p is a real number with p > 1, V (resp. V p ) stands for the set of 
probability measures on M (resp. with finite moment of order p) and dx for the Lebesgue 
measure on M. 

The Wasserstein distance of order p, valued in 1U {+oo}, is defined on V x V by 



where 7r runs over the set of probability measures on M? with marginals /z and \i. It is 
equivalently defined by 



where the infimum is taken over all random variables and X/x on the probability space 
(]0, 1[, dw) with respective laws /i and fi. It takes finite values on V p x V p and indeed 
defines a distance on V p . 

For complete references about the Wasserstein distances and related topics the reader 
can refer to [(]]. We only mention that both infima in ^ and © are achieved, and for 
the second definition we shall precise some random variables that achieve the infimum. 
For this purpose we introduce the notion of generalized inverse: 



W p (u x (t,.),u x (t,.)) = W p (u° x ,u° x ). 
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Definition 2.1 Let v belong to U . Then its generalized inverse is the function v 1 defined 
on ]0, 1[ by 

v~~ (w) = inf{x G M; v(x) > w}. 
Then t? _1 is a nondecreasing random variable on (]0, 1[, dw) by definition, with law v x 



since 



1 r l 



f(v 1 (w)) dw = J f'(s)l {s < v -i (w)} dsj dw 

f'(s)l{ v ( s )< w } dw) ds 



f'(s)(l-v(s)) ds 
f(s) dv x (s) 

for all / in C^(R). In particular its repartition function is v. 

Moreover this generalized inverse achieves the infimum in Q: 

Proposition 2.2 Let v and v in IA. Then we have (with possibly infinite values) 

'i - i/p 



W p (v x ,v x ) 



\v' 1 (w)-v- 1 (w)\ p dw\ 



for all p> 1. In particular for p = 1 we also have 

Wi(v x ,V x ) = \\V - 5||H(R). 

Proof. The general result is proved in The result specific to the case p — 1 follows 
by introducing, for a given v G U, the map defined on Mx]0, 1[ by 



1 if v(x) > w 
if v(x) < w, 



jv(x, w) = 
for which we have 

\v _1 — v^Kw) = / \jv — jv\(x, w) dx 
Jr 

for almost every w G ]0, 1[, and 

\jv — jv\(x,w) dw = \v — v\(x) 



for almost every i6M. Integrating the first equality on w in ]0, 1[ and the second one on 
x in E, we deduce 



i 

i „~,-i 



v —v | (it?) dw = / \v — v\(x) dx. 



o 
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Given v G U, its generalized inverse w _1 is actually the a.e. unique nondecreasing 
random variable on (]0, 1[, dw) with law v x . Given any other random variable X on 
(]0, l[,dw) with law v x , v~ l is called the (a.e. unique) nondecreasing rearrangement of X 
(see [H]). 

We conclude this section recalling a result relative to the convergence of probability 
measures. A sequence (fj, n ) of probability measures on R is said to converge weakly toward 

a probability measure fi if, as n goes to +oo, / ^pd\x n tends to if dp for all bounded 

Jr Jr 

continuous real-valued functions tp on R (or equivalently for all C°° functions tp with 
compact support, that is, if fi n converges to \x in the distribution sense). Given p > 1 this 
convergence is metrized on V p by the distance W p as shown by the following proposition 
(see pj): 

Proposition 2.3 Let p > 1, (fi n ) a sequence of probability measures in V p and /j, e V. 
Then the following statements are equivalent: 

i) (Wp(n n , //)) converges to 0; 

ii) (/i„) converges weakly to fi and sup / \x\ p dp, n (x) tends toO as R goes to infinity. 

n J\x\>R 

In this proposition we do not a priori assume that \i belongs to V p , but it can be noted 
that this property is actually induced by any of both hypotheses i) and ii). 
For measures in V we have the weaker result: 

Proposition 2.4 Let p > 1, (/i n ) and (u n ) two sequences in V converging weakly to fj, 
and v inV respectively. Then (with possibly infinite values) 

Wpin, v) < liminf W p (fi n , v n ). 

3 The case of classical solutions: Theorem H3 and corol- 
lary 



3.1 Proof of Theorem fTTTl 

We consider two classical solutions u and u to (P) such that u(t, .) and u(t, .) belong 
to U and are increasing for all t > 0, and we shall prove that 



W p (u x (t,.),u x (t,.)) = W p (u° x ,u° x ) 
as a consequence of Proposition 12.21 
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The map u° is increasing from to 1, so has a (true) inverse X(0, .) defined on ]0, 1[ 

by 

u°(X(0,w)) = w. 

Then, given w e]0, 1[, we consider a characteristic curve 1 1— > X(t,w) solution of 

(6) X t (t,w) = f'(u{t,X{t,w))) 

for t > 0, and taking value X(0, w) at t — 0. Since / is C 1 and u is bounded there exists a 
(non necessarily unique) solution X(.,w) to (JBJ by Peano Theorem (see [H] for instance); 
moreover by a classical computation from (QJ it is known to satisfy 

(7) u(t,X(t,w)) = w 
for allt > 0, from which it follows that 

X t (t,w) (= f(u(t,X(t,w))j)=f'(w) 

and hence 

(8) X(t,w) = X(0,w) + tf\w). 

In particular there exists a unique solution X(., w) to (JBJ). Now given t > 0, X(t, .) is the 
(true) inverse of the increasing function u(t, .) (by ((ZJ)), and Proposition 12.21 writes 



i/p 



W p (u x (t,.),u x (t,.)) = (jf |X(t,w)-X(t,ti;)|Pd«;) 
But from (JHJ) we obtain 

(9) X(*,u;) - X(t,w) = X(0,to) - X(0,w). 

This result ensures in particular that W p (u x (t, .),u x (t, .)) remains constant in time, may 
its initial value be finite or not; note however that (JSJ) is actually much stronger that 
Theorem O 



3.2 Proof of Corollary fTTTTI 

We assume that / is a C 2 convex function on K, and u° is a C 1 increasing initial profile 
in Up. 

First of all we note that the associated entropy solutions u is a classical solution in 
view of Theorem 11.21 this result is proved in [H] for instance, and its proof also ensures 
that u(t, .) is increasing for all t > 0. 
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Then we check that the moment property is preserved by the conservation law, that 
is, that u(t, .) also belongs to U p for any t > 0. Indeed, given t > 0, we have by the change 
of variable w = [u(t, -)]{x): 



which is finite since 



\X{t,w)\ p dw 
\X(0,w) + tf(w)\ p dw 



\X(0,w)rdw+F\\f\\loo m[) 



\X(0,w)\ p dw = / \x\ p u° x (x)dx 
Jr 

is finite by assumption. This ends the proof of Corollary II .111 



□ 



4 Time discretization of the conservation law 

In the previous section we have seen that the classical solutions are obtained through 
the method of characteristics, that we now summarize in our case: given an initial profile 
u° in U such that the corresponding solution u is C 1 and increasing in x for alH > 0, let 
X(0, .) be its inverse, defined by 

u°(X(0,w)) = w 
for all w e]0, 1[. Let then X(0,w) evolve into 
(10) X(t,w) = X(0,w)+tf'(w) 

for alH > and w 6 ]0, 1[ (see JEJ). The solution u(t, .) is then the inverse of the increasing 
map X(t, .), that is, is the unique solution of 

X(£, u(t, x)) = x. 



In the general case, defining X(0, .) in some similar way, there is no hope for the 
function X(t, .) defined by (fTTH) to be increasing for t > 0; inverting it would thus lead to 
a multivalued function, and no more to the entropy solution of the conservation law, as 
in the particular case discussed above. 

However, averaging (or "collapsing") this multivalued function into a single-valued 
function, Y. Brenier showed in how to build an approximate solution to the conservation 
law. 

We now precisely describe this so-called Transport-Collapse method in our case. 
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4.1 Definition and W p contraction property of the discretized so- 
lution 

Let u° G U be some fixed initial profile, with generalized inverse X(0, .) given as in 
Definition 12.11 by 

X(0, w) = inf{x G R; u°(x) > w} 

for all w g]0, 1[. X(0, .) can be seen as a random variable on the probability space ]0, 1[ 
equipped with the Lebesgue measure dw; its law is u^, ; sis pointed out after Definition 12 .11 
We let then X(0, .) evolve according to the method of characteristics, denoting 

X(h,w) = X(0,w) + hf(w) 

for all h > and almost every w g]0, 1[. Again, given h > 0, X(h,.) can be seen 
as a random variable on ]0, 1[; let then ThU° be its repartition function, that is, the 
function belonging to U and defined at any i e las the Lebesgue measure of the set 
{w g]0,1[; X(h,w) < x}. It is given by 

T h u°(x) = / l{x(h,w)<x}(w) dw. 
Jo 

We summarize this construction in the following definition: 
Definition 4.1 Let v G U with generalized inverse X(0, .) defined on ]0, l[ by 

X(0, w) = inf{x G R; v(x) > w}. 
Then, given h > 0, and letting 

X(h,w) = X(0,w) + hf'(w) 
for almost every w G ]0, 1[ ; we define the U function T^v on R by 

T h v(x) = / l {X (h,w)<x}(w) dw. 
Jo 

In the case of Section El (see (JHJ)), it turns out that X(h, .) is the (true) inverse of ThU°, 
and (h,x) i— > T h u°(x) is exactly the entropy solution to equation (JTJ with initial datum 
u° in U. This does not hold anymore in the general case, but will allow us to build an 
approximate solution Sh,u° by iterating the operator Th- Let us first give two important 
properties of Th'. 

Proposition 4.2 Let h>0,T h defined as above and p > 1. Then 

i) ThV belongs to U p if so does v. 

ii) For any v and v inU we have (with possibly infinite values unless v and v G U p ) 



W p ([T h v} x ,[T h v] x ) < W p (v x ,v x ). 
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Proof. It is really similar to what has been done in Section 01 as for Corollary II .111 
i) T h v belongs to U as a repartition function of a random variable, and we have 

\x\ p d[T h v\ x {x) = [ \X(h,w)\ p dw 



o 



i 

\X(0,w) + hf'(w)\ p dw 



o 



< 2 p - x [ \X(0,w)\ p + \hf(w)\ p dw 
Jo 

\x\ p dv x (x) + h p \\f H^ooqq^q 



< 2 p - x 

which ensures that [T^vj^ has finite moment of order p if so does v x . 

ii) On one hand the generalized inverses X(0, .) and X(0, .) of v and v respectively 
satisfy 

(11) W p (v x ,v x ) = (J\x(0,w)-X(0,w)\ p dwy /P 

by Proposition 12.21 (with finite values if both v and v belong to U p , and possibly infinite 
otherwise). On the other hand X(h, .) and X(h, .) have respective law [T^x and ph{)] x , 
so 



i/p 



(12) W p ([T hU °} x ,[T h u } x ) < (jf \X(h,w)-X(h,w)\ p dn?) 

by definition of the Wasserstein distance. But 

X(h,w)- X(h,w) = X(0,w) - 1(0, w) 

for almost every w G]0, 1[ by definition, which concludes the argument by (fTTf and (fT2|) . 

Note again that (|T2l holds only as an inequality since X(h, .) and X(h, .) are not 
necessarily nondecreasing, which was the case in the example discussed in Sectional □ 

We now use the operator Th defined above to build an approximate solution ShU° to 
the conservation law (HJ): 

Definition 4.3 Let h be some positive number and v G U. For any t > decomposed as 
t = (N + s)h with N G N and < s < I, we let 

S h v(t,.) = (l-s)T»v(.) + sT h N+ h(-) 
where T®v = v and T^ +1 v = T h {T^v). 
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These iterations make sense because T h v G U if v G U, Shv(t, .) EU (resp. U p ) for any 
h,t>0 and v & U (resp. U p ). 

We now prove two contractions properties on these approximate solutions. We first 
have the L X (R) contraction property: 

Proposition 4.4 Let h be some fixed positive number and Sh defined as above. Then, 
for any v and s, t > we have 

\\Shv(t, .) - S h v(s, - I* ~ s \ ll/'IU°°([o,i])- 

Proof. As in [T] we first observe that \\T h V — V\\ L i(^ < ^||/'||l<x.(]o,i[) for any V G U, 
then let t = (M + fi)h and s = (N + u)h with M, N G N and < )i, v < 1, and, for 
instance assuming that M > N, prove the proposition by applying this first bound to 
V= S h v(kh,.) = T£v for k = N + 1,...,M- 1. □ 

Then we have the W p contraction property: 

Proposition 4.5 Let h be some fixed positive number and Sh defined as above. Then, 
given v and v in U, we have for any t > 0: 

W p ([S h v] x (t,.),[S h v] x (t,.)) < W p (v x ,v x ). 

Proof. It follows from Proposition 14.21 (about T h ) and to the convexity of the W p 
distance to the power p, in the sense that 

WP(afix + (1 - a)fi 2 , av x + (1 - a)v 2 ) < aW*(ji!, + (1 - a)W^ 2 , v 2 ) 

for all real number a G [0,1] and probability measures ^i^ 2l V\ and v 2 (see |H] for in- 
stance). □ 

We shall now recall the convergence of the scheme toward the entropy solution of the 
conservation law. 

4.2 Convergence of the scheme in the Lj oc (M.) sense 

In this section we consider the space C([0, +oo[, L} oc (WL)) equipped with the topology 
defined by the semi-norms 



q nm (f) = sup / |/(ar)|da; 

te[0,n] J -m 

for any integers n and m and / G C([0, +oo[, Ll oc (R)). Then we have 

Proposition 4.6 Let u° G U. Then, as h goes to 0, the function ShU° converges in 
C([0, +oo[, Lj oc {R)) to the entropy solution of (GJj with initial datum u°. 
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We briefly give the steps of the proof, which follows the one of Brenier in [T], adapted 
to functions of U instead of L 1 (R). 

We first prove that the family (ShU°)h is relatively compact in C([0, +00 [, Lj oc (R)) by 
means of Proposition 031 and Helly and Ascoli-Arzela Theorems. Then we check that the 
limit of any sequence of (ShU°)h converging in C([0, +00 [, L} oc {R)) is an entropy solution 
to the conservation law (|T]) with initial datum u°. By the uniqueness of this solution 
ensured by Theorem II .![ this concludes the proof of Proposition 14.61 

4.3 Convergence of the scheme in W p distance sense 

We first prove a uniform equiintegrability result on the approximate solutions: 

Proposition 4.7 Let Sh be defined as above, v G U p and T > 0. Then 



sup sup / \x\ p d[Shv] x (t, x) 

)<h<T0<t<T J\x\>R 



0<h<T 0<t<T J \ X \>R 

tends to as R goes to infinity. 

Proof. We again denote M = ||/'||i,°o(]o,i[), and first consider Th itself, writing 

x\ p d[T h v] x (x) = / + hf(w)\ p l{\ v -i ( w )+hf( w )\>R} dw 

\x\>R JO 

/ X 



< / (\x\ + hM) p l { \ x \ +hM >R} dvj 

J 

< (l + J^Yf \A P dv x {x) 



R — hM J J\ x \>R-hM 
for R > hM. From this computation we deduce by iteration 

A' 



/ \x\ p d[T»v] x (x) < f[ f 1 



— rnrr) f \x\ p dv x (x) 

R-jhMj J lxl>R . NhM l 



for R > NhM, with 



3=1 

Thus 



fr/ hM \ ( hM V / NhM \ 

Uy' + R—JhMj-V + R-NhM) ~ GXP \R — NhM ) ' 



J WHW < ex P (^) / 



\x\ p dv x {x) 



-TM 



for any and h such that Nh < T. 
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>From this we get for instance 

/ 9«TA/f \ r 

\x\ p dv x (x) 



[ \x\>d[S h vUt t x) < exp (J^^) j 



\x\>R-2TM 

for any t and h smaller than T. This concludes the argument since the last integral tends 
to as R goes to infinity. □ 

>From this we deduce the convergence of the scheme in W p distance sense: 

Proposition 4.8 Let u° G U p and u be the entropy solution to (QJj with initial datum u . 
Then, for any t>0, W p ([ShU°] x (t, .),u x (t, .)) converges to as h goes to 0. 



Proof. Given t > 0, ShU°(t, .) converges to u(t, .) in Lj oc (M.) as h goes to (by 
Proposition I4.6jl . so [Shu] x (t, .) converges to the probability measure u x (t, .), first in the 
distribution sense, then in the weak sense of probability measures, and finally in W p 
distance by Propositions 14.71 and 12.31 

Note in particular that u x (t, .) has finite moment of order p for any t > 0, that is, 
u(t, .) belongs to U p . □ 

5 The general case of entropy solutions: Theorem II. HI 
and corollaries 



5.1 Proof of Theorem 11.61 

We let p > 1 and consider two initial data u° and u° in U with associated entropy 
solutions u and u. 

Given t > 0, Proposition 14.61 yields again the convergence of [ShU°] x (t, .) to u x (t, .) in 
the weak sense of probability measures. Since this holds also for u°, we obtain 

W p (u x (t, .),u x (t, .)) < liminf W p ([S h u°] x (t, .), [S h u°} x (t, .)) 

h—>0 

by Proposition 12.41 But, for each h, 

W p {[S h u°] x (t, .), [S h u°} x (t, .)) < W p (u° x ,u° x ) 
by Proposition 14.51 so finally 

W p {u x {t,.),u x {t,.))< W p (u° x ,u° x ). 
This concludes the argument. 
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5.2 Proof of Corollary fT9l 



We recall that in the introduction we have defined a distance on each U p by letting 



and we now prove that, given p > 1 and u £U P , the entropy solution u to the conservation 
law JTJ) belongs to C([0, +oo[,U p ). 

We first note, in view of the proof of Proposition 14.81 that u(t, .) indeed belongs to U v 
for all t > 0. 

Then, given s > 0, we need to prove that d p (u(t, .), u(s, .)) (= Wp(Ma.(t, .), u x (s, .))) 
tends to as t goes to s. Indeed, on one hand u(t, .) tends to u(s, .) in Lj oc (R) by 
Theorem II .1\ so u x (t, .) tends to u x (s, .), first in the distribution sense, then in the weak 
sense of probability measures. 

On the other hand, given T > s, we now prove that sup / \x\ p du x (t, x) goes to 

0<t<T J\ X \>R 

as R goes to infinity. For this, given e > 0, let R such that 



by Proposition EZI Let then yj G C C °°(R) such that < ip < 1 and y>(a;) = if |x| < R. 
On one hand 



as h goes to since tp(x)\x\ p G C^°(M) and [ShU°] x (t, .) tends to .) in distribution 
sense. On the other hand 



d p (u,u) = W p (u x ,u x ) 



sup sup / 

0<h<T 0<t<T J \ X \>R 






for all < h, t < T. Hence at the limit 




for all t <T, from which it follows that 




which means that indeed 




From these two results we deduce the continuity result by Proposition 12.31 
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5.3 Proof of Corollary fTTTOl 

Given t > 0, ShU°(t,.) converges to u(t, .) in L* oc (R) by Proposition 14.61 so for all 
s, t, n > we have 

||m(£, .) - u(s, .)\\mi-n,n}) = J™ II<Sa«°(*> •) ~ S h u°(8, OIUnt-M)- 

But 

||Sfcli°(t, .) - S h U°(s, .)|Ui([-n,n]) < •) - S/»li°(s, OIU 1 ^) - I* ~ S l II/'IU^W 

for all h > by Proposition H3J so letting h go to we get 

IK*, •) - OIUhi-™.™]) - I* _ s l ll/'IU°°(M)- 

Since this holds for all n > 0, we obtain Corollary II .101 

6 Extension to viscous conservation laws 

In this section we let v be a positive number and consider the viscous conservation 

law 

(13) Ut + f(u) x = vu xx t > 0, x e R 

with initial datum u° G L°°(R). 

Assuming that / is a locally Lipschitz real-valued function on M, and calling solution 
a function u in L°°([0, +oo[xIR) such that ([TT^ holds in the sense of distributions, it is 
known that, given u° G L°°(R), there exists a unique solution u to (|T3|) . If moreover 
n° G U, then w(t, .) also belongs to U for all t > 0, and the W p contraction property 
stated in Theorem 11.61 in the inviscid case v = still holds: 

Theorem 6.1 Given a locally Lipschitz real-valued function f on R and two initial data 
u° and u° in U, let u and u be the associated solutions to (fl3| . Then, for any t > and 
p > 1, we have (with possibly infinite values) 

W p (u x (t,.),u x (t,.)) < W p (u° x ,u° x ). 

We briefly mention how this contraction property for the viscous conservation law 
allows to recover the same property for the inviscid equation, given in Theorem 11.61 
Given some initial datum u° in IA and v > 0, let indeed u v be the corresponding solution 
to the viscous equation (TTTfl) . Then it is known (see jH| for instance) that u u (t, .) converges 
in L} oc ($t) to the solution u(t, .) to the inviscid conservation law (JTJ) with initial datum 
u°. From this the argument already used in Section I57T1 (with ShU°(t, .) intead of u u (t, .)) 
enables to recover Theorem II .61 
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The proof of Theorem 16 . II follows the lines of Sections 0] and Eland makes use of a time- 
discretization of equation (fl~3|) based on the discretization of the inviscid conservation law 
previously discussed. More precisely, given a time step h > 0, we first map u° e U to 
ThU° as in section 4.1, and then let ThU° evolve along the heat equation on a time interval 
h, that is, map it to 

T h u° = K h * T h u° 
where K h is the heat kernel defined on R by 



4,h 



Then, defining an approximate solution ShU° by iterating the % operator as in Definition 
14. 3| we prove that Propositions 14 . 21 and 14 . 51 still hold for the new 7^ and Sh operators (Note 
that the convolution with the heat kernel is a contraction for the Wasserstein distance of 
any finite order). 

Proposition HiU only holds assuming that v is twice derivable with v" in L X (R): it more 
precisely reads 

\\S h v(t, .) - S h v(s, .)|| L i(R) < \t-s\ [||/'||l-(]o,i[) + IK'IUhk)]- 

As in Section 14 .2\ this enables to prove that, given u° in U, twice derivable with (u )" 
in L 1 (R), the family (ShU°)h converges in C([0, +oo[, ^^(R)) to the solution of (fTH| with 
initial datum u°. 

With this convergence result in hand we follow the lines of Section 15.11 to prove The- 
orem 16.11 in the case of twice derivable initial data, with L 1 second derivative, while the 
general case follows by a density argument. 
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